# Importing Data
GERMANY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Germany.xlsx")
library(readxl)
Belgium <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Belgium.xlsx")
View(Belgium)
# Importing Data
GERMANY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Germany.xlsx")
# Creating Time Series Data
GERMANY_ts <- ts(GERMANY, start=c(2004,1), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
GERMANY_ts
sum(is.na(GERMANY_ts))
library(forecast)
GERMANY_ts <- tsclean(GERMANY_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(GERMANY_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
GERMANY_ts_model <- auto.arima(GERMANY_ts)
GERMANY_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
GERMANY_ts_forecast <- forecast (GERMANY_ts_model, level=c(95), h=257)
plot(GERMANY_ts_forecast)
GERMANY_ts_forecast
write.table(GERMANY_ts_forecast, file="Germany_TSA.csv", sep=",")
# Importing Data
JAPAN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Japan.xlsx")
# Creating Time Series Data
JAPAN_ts <- ts(JAPAN, start=c(2001,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
JAPAN_ts
sum(is.na(JAPAN_ts))
library(forecast)
JAPAN_ts <- tsclean(JAPAN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(JAPAN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
JAPAN_ts_model <- auto.arima(JAPAN_ts)
JAPAN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
JAPAN_ts_forecast <- forecast (JAPAN_ts_model, level=c(95), h=258)
plot(JAPAN_ts_forecast)
JAPAN_ts_forecast
write.table(JAPAN_ts_forecast, file="Japan_TSA.csv", sep=",")
# Importing Data
CHINA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/China.xlsx")
# Creating Time Series Data
CHINA_ts <- ts(CHINA, start=c(2005,1), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
CHINA_ts
sum(is.na(CHINA_ts))
library(forecast)
CHINA_ts <- tsclean(CHINA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(CHINA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
CHINA_ts_model <- auto.arima(CHINA_ts)
CHINA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
CHINA_ts_forecast <- forecast (CHINA_ts_model, level=c(95), h=257)
plot(CHINA_ts_forecast)
CHINA_ts_forecast
write.table(CHINA_ts_forecast, file="China_TSA.csv", sep=",")
# Importing Data
INDONESIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Indonesia.xlsx")
# Creating Time Series Data
INDONESIA_ts <- ts(INDONESIA, start=c(2006,01), end=c(2019,01), frequency=12)
# Viewing and Checking the Created Time Series Data
INDONESIA_ts
sum(is.na(INDONESIA_ts))
library(forecast)
INDONESIA_ts <- tsclean(INDONESIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(INDONESIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
INDONESIA_ts_model <- auto.arima(INDONESIA_ts)
INDONESIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
INDONESIA_ts_forecast <- forecast (INDONESIA_ts_model, level=c(95), h=263)
plot(INDONESIA_ts_forecast)
INDONESIA_ts_forecast
write.table(INDONESIA_ts_forecast, file="Indonesia_TSA.csv", sep=",")
# Importing Data
ITALY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Italy.xlsx")
# Creating Time Series Data
ITALY_ts <- ts(ITALY, start=c(2004,01), end=c(2018,12), frequency=12)
# Viewing and Checking the Created Time Series Data
ITALY_ts
sum(is.na(ITALY_ts))
library(forecast)
ITALY_ts <- tsclean(ITALY_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(ITALY_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
ITALY_ts_model <- auto.arima(ITALY_ts)
ITALY_ts_model
# Importing Data
ITALY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Italy.xlsx")
# Checking the Imported Data
View(ITALY)
# Creating Time Series Data
ITALY_ts <- ts(ITALY, start=c(2004,01), end=c(2018,12), frequency=12)
# Viewing and Checking the Created Time Series Data
ITALY_ts
sum(is.na(ITALY_ts))
library(forecast)
ITALY_ts <- tsclean(ITALY_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(ITALY_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
ITALY_ts_model <- auto.arima(ITALY_ts)
ITALY_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
ITALY_ts_forecast <- forecast (ITALY_ts_model, level=c(95), h=264)
plot(ITALY_ts_forecast)
ITALY_ts_forecast
write.table(ITALY_ts_forecast, file="Italy_TSA.csv", sep=",")
# Importing Data
BELGIUM <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Belgium.xlsx")
# Creating Time Series Data
BELGIUM_ts <- ts(BELGIUM, start=c(2004,1), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
BELGIUM_ts
sum(is.na(BELGIUM_ts))
library(forecast)
BELGIUM_ts <- tsclean(BELGIUM_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(BELGIUM_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
BELGIUM_ts_model <- auto.arima(BELGIUM_ts)
BELGIUM_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
BELGIUM_ts_forecast <- forecast (BELGIUM_ts_model, level=c(95), h=258)
plot(BELGIUM_ts_forecast)
BELGIUM_ts_forecast
write.table(BELGIUM_ts_forecast, file="Belgium_TSA.csv", sep=",")
# Importing Data
CANADA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Canada.xlsx")
# Creating Time Series Data
CANADA_ts <- ts(CANADA, start=c(2007,1), end=c(2019,03), frequency=12)
# Viewing and Checking the Created Time Series Data
CANADA_ts
sum(is.na(CANADA_ts))
library(forecast)
CANADA_ts <- tsclean(CANADA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(CANADA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
CANADA_ts_model <- auto.arima(CANADA_ts)
CANADA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
CANADA_ts_forecast <- forecast (CANADA_ts_model, level=c(95), h=261)
plot(CANADA_ts_forecast)
CANADA_ts_forecast
write.table(CANADA_ts_forecast, file="Canada_TSA.csv", sep=",")
# Importing Data
FRANCE <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/France.xlsx")
# Creating Time Series Data
FRANCE_ts <- ts(FRANCE, start=c(2004,1), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
FRANCE_ts
sum(is.na(FRANCE_ts))
library(forecast)
FRANCE_ts <- tsclean(FRANCE_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(FRANCE_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
FRANCE_ts_model <- auto.arima(FRANCE_ts)
FRANCE_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
FRANCE_ts_forecast <- forecast (FRANCE_ts_model, level=c(95), h=257)
plot(FRANCE_ts_forecast)
FRANCE_ts_forecast
write.table(FRANCE_ts_forecast, file="France_TSA.csv", sep=",")
# Importing Data
SPAIN <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Spain.xlsx")
# Creating Time Series Data
SPAIN_ts <- ts(SPAIN, start=c(2004,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
SPAIN_ts
sum(is.na(SPAIN_ts))
library(forecast)
SPAIN_ts <- tsclean(SPAIN_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(SPAIN_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
SPAIN_ts_model <- auto.arima(SPAIN_ts)
SPAIN_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
SPAIN_ts_forecast <- forecast (SPAIN_ts_model, level=c(95), h=257)
plot(SPAIN_ts_forecast)
SPAIN_ts_forecast
write.table(SPAIN_ts_forecast, file="Spain_TSA.csv", sep=",")
# Importing Data
SLOVAKIA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Slovakia.xlsx")
# Creating Time Series Data
SLOVAKIA_ts <- ts(SLOVAKIA, start=c(2004,05), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
SLOVAKIA_ts
sum(is.na(SLOVAKIA_ts))
library(forecast)
SLOVAKIA_ts <- tsclean(SLOVAKIA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(SLOVAKIA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
SLOVAKIA_ts_model <- auto.arima(SLOVAKIA_ts)
SLOVAKIA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
SLOVAKIA_ts_forecast <- forecast (SLOVAKIA_ts_model, level=c(95), h=258)
plot(SLOVAKIA_ts_forecast)
SLOVAKIA_ts_forecast
write.table(SLOVAKIA_ts_forecast, file="Slovakia_TSA.csv", sep=",")
# Importing Data
NEITHERLAND <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Neitherland.xlsx")
# Creating Time Series Data
NEITHERLAND_ts <- ts(NEITHERLAND, start=c(2004,03), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
NEITHERLAND_ts
sum(is.na(NEITHERLAND_ts))
library(forecast)
NEITHERLAND_ts <- tsclean(NEITHERLAND_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(NEITHERLAND_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
NEITHERLAND_ts_model <- auto.arima(NEITHERLAND_ts)
NEITHERLAND_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
NEITHERLAND_ts_forecast <- forecast (NEITHERLAND_ts_model, level=c(95), h=258)
plot(NEITHERLAND_ts_forecast)
NEITHERLAND_ts_forecast
write.table(NEITHERLAND_ts_forecast, file="Neitherland_TSA.csv", sep=",")
# Importing Data
UK <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/UK.xlsx")
# Creating Time Series Data
UK_ts <- ts(UK, start=c(2004,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
UK_ts
sum(is.na(UK_ts))
library(forecast)
UK_ts <- tsclean(UK_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(UK_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
UK_ts_model <- auto.arima(UK_ts)
UK_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
UK_ts_forecast <- forecast (UK_ts_model, level=c(95), h=257)
plot(UK_ts_forecast)
UK_ts_forecast
write.table(UK_ts_forecast, file="UK_TSA.csv", sep=",")
# Importing Data
TURKEY <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Turkey.xlsx")
# Creating Time Series Data
TURKEY_ts <- ts(TURKEY, start=c(2002,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
TURKEY_ts
sum(is.na(TURKEY_ts))
library(forecast)
TURKEY_ts <- tsclean(TURKEY_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(TURKEY_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
TURKEY_ts_model <- auto.arima(TURKEY_ts)
TURKEY_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
TURKEY_ts_forecast <- forecast (TURKEY_ts_model, level=c(95), h=257)
plot(TURKEY_ts_forecast)
TURKEY_ts_forecast
write.table(TURKEY_ts_forecast, file="Turkey_TSA.csv", sep=",")
# Importing Data
THAILAND <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Thailand.xlsx")
# Creating Time Series Data
THAILAND_ts <- ts(THAILAND, start=c(2007,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
THAILAND_ts
sum(is.na(THAILAND_ts))
library(forecast)
THAILAND_ts <- tsclean(THAILAND_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(THAILAND_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
THAILAND_ts_model <- auto.arima(THAILAND_ts)
THAILAND_ts_model
THAILAND_ts_model <- auto.arima(THAILAND_ts, ic="aic", trace = TRUE)
THAILAND_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
THAILAND_ts_forecast <- forecast (THAILAND_ts_model, level=c(95), h=257)
plot(THAILAND_ts_forecast)
THAILAND_ts_forecast
write.table(THAILAND_ts_forecast, file="Thailand_TSA.csv", sep=",")
# Importing Data
KOREA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/Korea.xlsx")
# Creating Time Series Data
KOREA_ts <- ts(KOREA, start=c(2008,01), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
KOREA_ts
sum(is.na(KOREA_ts))
library(forecast)
KOREA_ts <- tsclean(KOREA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(KOREA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
KOREA_ts_model <- auto.arima(KOREA_ts)
KOREA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
KOREA_ts_forecast <- forecast (KOREA_ts_model, level=c(95), h=258)
plot(KOREA_ts_forecast)
KOREA_ts_forecast
write.table(KOREA_ts_forecast, file="Korea_TSA.csv", sep=",")
# Importing Data
USA <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/USA.xlsx")
# Creating Time Series Data
USA_ts <- ts(USA, start=c(2005,01), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
USA_ts
sum(is.na(USA_ts))
library(forecast)
USA_ts <- tsclean(USA_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(USA_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
USA_ts_model <- auto.arima(USA_ts)
USA_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
USA_ts_forecast <- forecast (USA_ts_model, level=c(95), h=257)
plot(USA_ts_forecast)
USA_ts_forecast
write.table(USA_ts_forecast, file="USA_TSA.csv", sep=",")
library(readxl)
France <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/France.xlsx")
View(France)
# Importing Data
FRANCE <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/France.xlsx")
# Creating Time Series Data
FRANCE_ts <- ts(FRANCE, start=c(2004,1), end=c(2019,07), frequency=12)
# Viewing and Checking the Created Time Series Data
FRANCE_ts
sum(is.na(FRANCE_ts))
library(forecast)
FRANCE_ts <- tsclean(FRANCE_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(FRANCE_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
FRANCE_ts_model <- auto.arima(FRANCE_ts)
FRANCE_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
FRANCE_ts_forecast <- forecast (FRANCE_ts_model, level=c(95), h=257)
plot(FRANCE_ts_forecast)
FRANCE_ts_forecast
write.table(FRANCE_ts_forecast, file="France_TSA.csv", sep=",")
# Importing Data
FRANCE <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/France.xlsx")
# Creating Time Series Data
FRANCE_ts <- ts(FRANCE, start=c(2004,1), end=c(2019,06), frequency=12)
# Viewing and Checking the Created Time Series Data
FRANCE_ts
sum(is.na(FRANCE_ts))
library(forecast)
FRANCE_ts <- tsclean(FRANCE_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(FRANCE_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
FRANCE_ts_model <- auto.arima(FRANCE_ts)
FRANCE_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
FRANCE_ts_forecast <- forecast (FRANCE_ts_model, level=c(95), h=258)
plot(FRANCE_ts_forecast)
FRANCE_ts_forecast
# Importing Data
FRANCE <- read_excel("C:/Users/Biniam/Desktop/Documents/Academic/Thesis/Analysis Folder/Excel Files/Lebanon/France.xlsx")
# Creating Time Series Data
FRANCE_ts <- ts(FRANCE, start=c(2004,1), end=c(2019,04), frequency=12)
# Viewing and Checking the Created Time Series Data
FRANCE_ts
sum(is.na(FRANCE_ts))
library(forecast)
FRANCE_ts <- tsclean(FRANCE_ts)
# Step – 1 of the Box-Jenkins Methodology (Identification: Plotting the Time Series Data)
plot(FRANCE_ts)
# Step-2 of the Box-Jenkins Methodology (Estimating the appropriate model)
FRANCE_ts_model <- auto.arima(FRANCE_ts)
FRANCE_ts_model
# Forecasting
options(max.print=1000000)
library(forecast)
FRANCE_ts_forecast <- forecast (FRANCE_ts_model, level=c(95), h=260)
plot(FRANCE_ts_forecast)
FRANCE_ts_forecast
write.table(FRANCE_ts_forecast, file="France_TSA.csv", sep=",")
